A new technique for the sp 2 /sp 3 characterisation of carbon materials 
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We present a technique to determine the sp 3 /sp 2 ratio of carbon materials which is based on 
the electron energy-loss spectroscopy and which uses the theoretical spectrum of graphite obtained 
from ab initio electronic structure calculations. The method relies on the separation of the tt* and 
a* components of the carbon K-edge of graphite. This tt* spectrum is adopted and assumed to be 
transferable to other carbon systems given an appropriate parametrisation of the broadening. The 
method is applied on a series of Monte-Carlo generated amorphous carbon structures and is shown 
to be stable over a wide range of the energy windows for which spectral integration is performed. 
The results are found to be in good agreement with the sp 3 fraction obtained from a microscopic 
scheme which uses the tt— orbital axis vector (POAV1) analysis. The technique was also applied on 
a series of experimental spectra of amorphous carbon and was found to be in good agreement with 
the results obtained from a functional fitting approach. 



One of the parameters that determine the hardness of 
amorphous carbon is the fraction of sp 3 bonded carbon 
atoms. This fraction can be determined by making use 
of the X-ray photoelectron spectroscopy (XPS) Q, the 
plasmon energy technique 0, Q] , the near-edge X-ray ab- 
sorption spectroscopy (NEXAFS) Q and the high-energ 
electron energy-loss spectroscopy (HEELS) diflST 

The high-energy EEL spectrum of an amorphous car- 
bon is made of a tt* part whose onset is at about 284 eV 
and a a* component whose onset is about 5-7 eV from 
that of the tt* component. The HEELS sp 3 /sp 2 charac- 
terisation technique is based on the ability to isolate the 
7r* features and subsequently its cross-section from the 
entire spectrum. All the isolation methods that are of- 
ten used are limited by the assumption that the tt* and 
the a* bands are well separated in energy such that their 
cross-sections can be independently deduced by spectral 
integrations over fixed energy windows. But it is well 
known that plural inelastic scattering such as core-loss ex- 
citation plus additional plasmon losses modifies the edge 
intensities by spreading oscillator strengths to higher en- 
ergies The tt* band can spread to energies as high 
as 30 eV beyond the edge onset thereby overlapping with 
the a* band. Density functional theory (DFT) calcula- 
tions on graphite and amorphous carbon, which will be 
presented elsewhere , reveal that the tt* band is excep- 
tionally broad, extending even to energies beyond 30 eV. 
Therefore, a sp 3 /sp 2 ratio characterisation that ignores 



these finite tt* intensities (of the HEELS or NEXAFS) 
at higher energies is not entirely consistent. In this let- 
ter we present a new technique for the characterisation of 
the sp 3 /sp 2 ratio in (amorphous) carbon systems which is 
based on the tt* /a* separation of the calculated energy- 
loss near-edge structure (ELNES) of graphite. We show 
that this consistent treatment of the 7r* and the a* spec- 
tra can lead to sp 3 proportions that are independent of 
the choice of the size of the integration energy window in 
sharp contrast with other techniques 

DFT calculations are performed using the all- 
electron Full-Potential-Linearised- Augmented-Plane- 
Wave (FLAPW) code WIEN2k |lfj. The exchange and 
correlation energy is treated using the local density 
approximation Muffin-tin radii (Rmt) are fixed 

at 1.3 atomic units (a.u.) for graphite and 1.25 a.u. 
for amorphous carbon, while RMTxK mal values of 5.5 
and 5.0 are used for graphite and amorphous carbon, 
respectively; K max being the plane wave cut-off. Up to 
1000 k-points are used to sample the full Brillouin zone 
of graphite (72 in the irreducible Brillouin zone (IBZ)) 
while 20 (4 in the IBZ) are used for the amorphous 
carbon. The core-hole effect is introduced in the ELNES 
of graphite via the so-called excited-core method (in a 
supercell approach). The tt* component of the ELNES 
of graphite is obtained using the expressions of ref. 
for the differential cross-section which was recently 
implemented on hexagonal boron nitride [l^ and on 
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nanotubes (l^ . 

In order to obtain insight into the local structure of 
amorphous carbon, we generate amorphous carbon sys- 
tems by means of classical Monte Carlo (MC) simula- 
tions based on the Tersoff empirical potential for 
carbon Q| . We generated amorphous carbon structures 
at densities of 2.0, 2.6, 3.0 and 3.6 g/cm 3 by quenching 
hot gaseous carbon from 8000K down to 300K. This was 
done via a continuous space, periodic boundary condition 
(PBC), constant volume, constant temperature and con- 
stant number of atoms classical Metropolis Monte Carlo 
procedure on cubic cells made of 64 atoms. The ELNES 
spectra (without core-hole) were computed for the four 
generated carbon systems. The average (over all atoms) 
ELNES spectra of each of the systems were calculated. 
On Figure 1 we show the ELNES spectra of the four gen- 
erated carbon systems. 
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Figure 1: Fotal ELNES (in the sudden approximation) of four 
MC-generated carbon structures of varying densities of 2.00- 
3.60 g/cm 3 . 



The simulated structures allow us to look at the lo- 
cal environment of every atom and to determine the sp 3 
fraction from a microscopic point of view. For curved 
carbon systems like the fullerenes and the nanotubes, 
the effect of curvature on the hybridisation state of car- 
bon atoms has been successfully described by u sing the 
7r-orbital axis vector (POAV1) analysis \v\ 111 lia l^of . 



The application of this technique to all three-coordinated 
atoms of our MC-generated structures provides a more 
precise approach to characterise the sp 3 fraction because 
it accounts for the non-planarity of the three-coordinated 
atoms. However, not all four-coordinated atoms are hav- 
ing the ideal sp 3 local bonding structure as follows from 
the local geometry and our calculations of the core-loss 
ELNES spectra of the generated structures. We can eas- 
ily account for the distortion-induced sp 2 character in 
such four-coordinated atoms by analysing the tt* com- 
ponents for every non ideal tetrahedral four-coordinated 
atom as explained in ref. 9]. The sp 3 fraction obtained 
by this approach is consistently higher than that obtained 
by using the coordination number method for low and in- 
termediate densities (see Table HJ. For high densities the 
two methods yield the same results. 

The technique of the sp 3 /sp 2 characterisation pre- 
sented here is based on the assumption that the ELNES 
of any carbon atom can be decomposed into two inde- 
pendent parts. One of these ELNES components (the tt* 
component) is sought by performing first-principles cal- 
culations on graphite for which the two components are 
well defined. To some extent, the method is compara- 
ble to the experimental work of ref. |2lJ. The sp 3 /sp 2 
quantification procedure is described as follows 

• We identify the energy position E w of the main tt* 
peak from the spectrum of the uncharacterised sam- 
ple and rescale the intensity to have a value of one 
at this peak. We denote this new spectrum as 1(E) 
where E is the energy loss. 

• We consider a normalised smoothening func- 
tion g(a, E) ( which could be a Gaussian or a 
Lorentzian) of full width at half maximum a. 

• We define an energy scale shift parameter s by 
which the unbroadened calculated graphite tt* spec- 
trum should be shifted to align the maximum of its 
main tt* peak to that of the unknown sample. After 
the spectral shift the tt* spectrum is scaled to one 
at E w . 
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• Using this new shifted graphite it* spectrum which 
is convolved with g(s,a,E), denoted as Tl(s,a, E), 
we define the function 

F(s,<r)= / (I(E)-U(s,a,E)) 2 dE (1) 

JE min 

where E m j„ is the edge-onset. In Eq.Q) it is im- 
plied that the spectrum of the amorphous carbon 
is entirely of tt* character up to the main tt* peak. 

• F is minimised for s and a. The minimal values 
of s and a are used for U(s, a, E) to deduce the 7r* 
spectrum which best fits the lower energy part of 
the unknown spectrum. This shifted and convolved 
it* spectrum is assumed to be that of the unknown 
amorphous carbon system. 

• An energy window AE whose lower limit is 
fixed at the edge onset is chosen and the in- 
tegrals I W (AE) = J^"" +AE U(s,a,E)dE and 
iTotai(AE) = fg m ' n+AE I(E)dE are evaluated. 



• The ratios, for various AE, R(AE) = 
I n (AE) I 'lTotai{AE) are calculated. 

• The ratio Ro(AE) for a well characterised sample, 
whose sp 2 fraction is x$ is known, is also calculated 
and the sp 2 fraction in the unknown amorphous 
carbon sample is deduced as 



R(AE) 

and the sp 3 fraction is 1 — x. 



(2) 



To test and validate this method, we applied this tech- 
nique to the ELNES spectra of the four generated amor- 
phous carbon systems described above. The ratio R(AE) 
for each system is plotted as a function of AE on Figure 2. 
From this figure it is seen that for each energy window 
R decreases monotonically as the density increases. It 
is also seen that for the wide range of AE considered 
(15 eV < AE < 55 eV) all the curves showing the varia- 
tion of R with AE seem to be parallel to each other. This 
suggests that if one of the spectra was a reference spec- 
trum then the ratio of sp 2 carbon in each spectrum with 



respect to this reference will be fairly constant over this 
wide energy range. This demonstrates that the method 
should yield sp 2 or sp 3 fractions that are independent 
of the choice of the energy window in contrast to other 
techniques. In those techniques the sp 3 fraction depends 
crucially on both the width and position of the energy 
window. 
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Figure 2: The ratio R=Area(7r*)/Area(7r*+cr*) for the seven 
generated carbon structures. Remark that all the lines are 
approximately parallel to each other. 



To determine the sp 3 fraction, we take one of the 
Monte-Carlo generated amorphous carbon structures as 
a reference. Because the sp 3 fraction determined by the 
POAV1 analysis coincides with that obtained by the near- 
est neighbour counting method for high densities (see Ta- 
bleQJ, the structure having a density of 3.6 g/cm 3 with an 
sp 3 fraction of 90% is considered as the reference. The re- 
sulting fractions for the other three systems are shown on 
Figure 3. The sp 3 fraction is clearly seen to be indepen- 
dent of the size of the energy window. The largest vari- 
ation recorded did not exceed 6%. The values obtained, 
reproduced in Tabled were in very good agreement with 
those obtained from a POAVl-based analysis. These re- 
sults are given as intervals because of the standard de- 
viations on the POAVl-corrected sp 3 fractions recorded 
on the entire systems of the three-coordinated atoms. In 
the table we also report the proportion of 4-coordinated 
atoms in each of the systems. 
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Figure 3: The sp 3 percentage for the four generated carbon 
structures. The p = 2>.6g/cm z sample made of 10.3% sp 2 
atoms (according to the POAV1 approximation) is used as 
the reference system. 
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Figure 4: The sp 3 percentage for the five amorphous carbon 
samples. 

Table|n]compares the sp 3 percentages of the five amor- 
phous carbon samples with results obtained using the 
functional fitting method. 



density 


new method 


coordination number 


POAV1 




Sample 


new method 


Fitting technique 


2.00 


32-38 


16 


28-32 




I 


58-68 


67-68 


2.60 


39-40 


20 


35-37 




II 


60-63 


61 


3.00 


56-58 


55 


54-56 




III 


56-57 


56 


3.60 


90 


90 


90 




IV 


54-55 


51-53 










V 


52-55 


53 



Table I: The sp percentages of the four generated amorphous 
carbon systems. The 3.6 g/cm 3 system is considered as ref- 
erence and its sp 2 ratio (0.103) according to the POAV1 is 
adopted. 



Table II: The sp percentages of the five amorphous carbon 
samples. The 3.6 g/cm 3 system is considered as reference 
and its sp 2 ratio (0.103) according to the POAV1 is adopted. 
Apart from sample I, all samples show excellently stable re- 
sults within the 30 eV window. 



The use of this technique is not restricted to computer 
generated samples (in contrast to the POAV1 analysis) 
but can also be applied to experimental spectra. On Fig- 
ure 4 and in Table |n] we show the sp 3 fractions of five 
spectra obtained from ref. ,22]. In this table we also 
listed the values obtained in |22| using a functional fit- 
ting method with three Gaussians. Remark the good 
agreement with the functional fitting technique. Sample 
I showed the strongest dependence on AE, but for the 
large energy window range of 15 eV - 45 eV the variation 
did not exceed 10% for this sample. 



In summary, We have developed a method to quantify 
the sp 3 fraction of carbon materials based on the the- 
oretical separation of the it* and a* components of the 
carbon Is electron energy-loss spectra. The 7r* spectrum 
of graphite is adopted and assumed to be transferable to 
other carbon systems given an appropriate parametrisa- 
tion of the broadening. This has enabled the isolation 
of the 7T* spectrum of the carbon systems and a subse- 
quent characterisation of their sp 3 contents. We have 
applied the method to a series of Monte-Carlo generated 
amorphous carbon structures and shown that the method 
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gives stable results over a wide range of the energy win- 
dow used for spectral integration and that the sp 3 frac- 
tions obtained using this method agree with values re- 
sulting from an exhaustive analysis of the distorsion- 
induced hybridisation of every carbon atom by means 
of the 7T— orbital axis vector (POAV1) technique. The 
method has also been applied to five experimental amor- 
phous carbon spectra and found to be in good agreement 
with the results obtained from a more traditional func- 



tional fitting technique. 
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